GoSam: A program for automated one-loop 
Calculations 



G. Cullen 

Deutsches Elektronen-Synchrotron DESY, Zeuthen, Germany 

N. Greiner, G. Heinricl:^^? T. Reiter 

Max-Planck-Institute for Physics, Munich, Germany 

G. Luisoni 

Institute for Particle Physics Phenomenology, University of Durham, UK 
P. Mastrolia 

Ma:x;-Planck Institute for Physics, Munich, Germany; 
Dipartimento di Fisica, Universita di Padova, Italy 

G. Ossola 

New York City College of Technology, City University of New York 

F. Tramontano 

CERN, Geneva, Switzerland 

Abstract. The program package GoSam is presented which aims at the automated calculation 
of one-loop amplitudes for multi-particle processes. The amplitudes are generated in terms of 
Feynman diagrams and can be reduced using either D-dimensional integrand-level decomposition 
or tensor reduction, or a combination of both. GoSam can be used to calculate one-loop 
corrections to both QCD and electroweak theory, and model files for theories Beyond the 
Standard Model can be linked as well. A standard interface to programs calculating real 
radiation is also included. The flexibility of the program is demonstrated by various examples. 



1. Introduction 

Precise tiieory predictions play an important role in the analysis and the interpretation of collider 
physics data, in the search for new particles as well as to constrain model parameters at a 
later stage. Therefore it is desirable to have predictions at next-to-leading order (NLO) in 
perturbation theory as a standard at the LHC. However, this can only be achieved if a level of 
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automation for NLO predictions is reached which comes close to the one we have for leading 
order tools. 

The need for an automation of NLO calculations has been noticed some time ago 
and lead to public programs like FeynArts [IJ and QGraf [2j for diagram generation and 
FormCalc/Loop Tools [3j and GRACE [2] for the automated calculation of NLO corrections. 
However, calculations of one-loop amplitudes with more than four external legs were still tedious 
case- by-case calculations. Only very recently, conceptual and technical advances in multi- 
leg one-loop calculations opened the door to the possibility of an automated generation and 
evaluation of multi-leg one-loop amplitudes. As a consequence, public NLO tools containing 
a collection of hard-coded individual processes, like e.g. MCFM 1^, VBFNLO [3 El E], 
MC@NLO [IDIIII], POWHEG-Box [IJitlSj, POWHEL [HIISIIIS], are now being complemented 
by flexible automated tools [T71 [TSl EHl [201 EI] which have as their aim that basically any process 
which may turn out to be important for the comparison of LHC findings to theory can be 
evaluated at NLO accuracy. 

In this talk, we present the program package GoSam[T7] which allows the automated 
calculation of one- loop amplitudes for multi-particle processes. Amplitudes are expressed in 
terms of Feynman diagrams, where the integrand is generated analytically using QGRAF [2], 
FORM [22j, spinney [23j and haggles [23]. The individual program tasks are steered via python 
scripts, while the user only needs to edit an "input card" to specify the details of the process to 
be calculated, and launch the generation of the source code and its compilation, without having 
to worry about internal details of the code. 

The program offers the option to use different reduction techniques: either the unitarity-based 
integrand reduction [251 ES] as implemented in Samurai [27J or traditional tensor reduction as 
implemented in golem95C |281 129] interfaced through tensorial reconstruction at the integrand 
level |30j . or a combination of both. It can be used to calculate one-loop corrections within 
both QCD and electroweak theory. Beyond the Standard Model theories can be interfaced using 
FeynRules [31] or LanHEP [32] . The Binoth Les Houches interface [33] to programs providing the 
real radiation contributions is also included. 

2. Theoretical framev^fork 

2. 1 . Generation of the diagrams 

For the diagram generation we use the program QGRAF [2j and supplement it by adding another 
level of analysing and filtering the diagrams, written in Python. This allows for example to drop 
diagrams whose colour factor turns out to be zero, or to determine the signs for diagrams with 
Majorana fermions. 

Within our framework, QGRAF generates three sets of output files: an expression for each 
diagram to be processed with FORM [22j, Python code to draw all diagrams, and Python code to 
compute the properties of each diagram. The information about the model is either read from 
the built-in Standard Model file or is generated from a user defined LanHEP [32] or Universal 
FeynRules Output (UFO) ^ file. 

The program also produces a LTgX file which contains graphical representations of all 
diagrams as well as a summary of the helicity and colour basis used. 

2.2. Code generation 

The amplitude is generated in terms of algebraic expressions based on Feynman diagrams and 
then processed with a FORM program, using spinney [23J for the spinor algebra. In GoSam we 
have implemented the 't Hooft-Veltman scheme (HV) and dimensional reduction (DRED). In 
both schemes all external vectors (momenta and polarisation vectors) are kept in four dimensions, 
while internal vectors are kept in the n-dimensional vector space (n = 4 — 2e). We adopt the 
conventions used in [2^, where k denotes the four-dimensional projection of an n-dimensional 



vector k. The (n— 4)-dimensional orthogonal projection is denoted by k. For the loop momentum 
q we introduce the symbol /x^ = —(f, such that 

q =q +q =q - H . (1) 

To prepare the numerator functions of the one- loop diagrams for their numerical evaluation, 
we separate the symbol //^ and dot products involving the momentum q from all other factors. 
All subexpressions which do not depend on either q or ^u^ are substituted by abbreviations, 
which are evaluated only once per phase space point. Each of the two parts is then processed 
by haggles [24j, which generates optimised Fortraii95 code for the numerical evaluation. For 
each diagram we generate an interface to Samurai [27j, goleiii95C [29] and/or PJFRY [34 1. 

Our standard choice for the reduction is to use Samurai [27], which provides a very fast and 
stable reduction in a large part of the phase space. Furthermore, Samurai reports to the client 
code if the quality of the reconstruction of the numerator suffices the numerical requirements. In 
GoSam we use this information to trigger an alternative reduction with either golem95C [29] or 
PJFRY [34J whenever these reconstruction tests fail. This combination of on-shell techniques and 
traditional tensor reduction is achieved using tensorial reconstruction at the integrand level [30] . 
The tensorial reconstruction not only can cure numerical instabilities, but in some cases also can 
reduce the computational cost of the reduction. Since the reconstructed numerator is typically 
of a form where kinematics and loop momentum dependence are already separated, the use of 
a reconstructed numerator tends to be faster than the original procedure, in particular in cases 
with a large number of legs and low rank. 

2.3. Rational terms 

Terms containing the symbols fi^ or e in the numerator of the integrands can lead to a so-called 
i?2 term [35J, which contributes to the rational part of the amplitude. We generate the R2 part 
along with all other contributions without the need to seperate the different parts. In addition, 
we provide an implicit and an explicit construction of the R2 terms, using the fact that there 
are two ways of splitting the numerator function: 

M{q,^?,e) = AAo(g,//2)+eAAi(g,/i2) + e2_^2(g,^') (2) 
or, alternatively, 

^^{q,^i\e) = M{q)+M{q,i2^e). (3) 

The implicit construction uses the splitting of Eq. ^ and treats all three numerator functions J\fi 
on equal grounds. Each of the three terms is reduced seperately in a numerical reduction and 
the Laurent series of the three results are added up taking into account the powers of e. 

The explicit construction of R2 is based on the assumption that each term in J\f in Eq. ([3]) 
contains at least one power of //^ or e. The expressions for those integrals are relatively simple 
and known explicitly. Hence, the part of the amplitude which originates from TV is computed 
analytically whereas the purely four-dimensional part Af is passed to the numerical reduction. 

In the program, possible options for R2 are r2=iiiiplicit, explicit, off and only. Using 
r2=only discards everything but the R2 term and puts GoSam in the position of providing R2 
terms to supplement other codes which work entirely in four dimensions. 

2.4. Conventions 

To be specific, we consider the case where the user wants to compute QCD corrections. In the 
case of electroweak corrections, the analogous conventions apply except that the strong coupling 
gs is replaced by e. In the QCD case, the tree-level matrix element squared can be written as 

\M\i,, = AlAo = {9sf'-ao, (4) 



where 6 = is also possible. The matrix element at one-loop level, i.e. the interference term 
between tree-level and one-loop, can be written as 



lA^liioop = Aa, + AIA, = 2.^{AlA,) = {9s?' h + ^ + ^ + 0{e) 



(5) 



A call to the subroutine Samplitude returns an array consisting of the four numbers 
(oq, Co, c_i, c_2) in this order. The average over initial state colours and helicities is included in 
the default setup. Renormalisation is included depending on the options chosen by the user, for 
a more detailed description we refer to [T7]. In cases where the process is loop induced, i.e. the 
tree level amplitude is absent, the program returns the values for 

asifi) (47r)^ 
271 T{l-e] 

has been extracted. 



3. Installation and Usage 

3.1. Installation 

The user can download the code GoSam either as a tar-ball or from the subversion repository 
at 

http : //projects .hepforge . org/gosam/ . 

The build process and installation of GoSam is controlled by Python Distutils, while the 
build process for the libraries Samurai and golem95C is controlled by Autotools. To install 
GoSam, the user needs to run 

python setup. py install — prefix MYPATH 

If MYPATH is different from the system default, the environment variables PATH, LD .LIBRARY _PATH 
and PYTHONPATH might have to be set accordingly. For more details we direct the user to the 
GoSam reference manual coming with the code. 

The program is designed to run in any modern Linux/Unix environment; we expect that 
Python (> 2.6), Java (> 1.5) and Make are installed on the system. Furthermore, a Fortran 95 
compiler is required in order to compile the generated code. 

On top of a standard Linux environment, the programs FORM |22j, version > 3.3, and 
QGRAF [2] need to be installed on the system. Further, at least one of the libraries Sa- 
murai [27] and golem95C [29] needs to be present at compile time of the generated code. 
For the user's convenience we have prepared a package containing Samurai and golem95C 
together with the integral libraries OneLOop |36j. QCDLoop [37j and FF[38j. The package 
gosam-contrib-1 . . tar .gz containing all these libraries is available for download from 
[http : //pro j ects . hepforge . org/gosam/ 

3.2. Usage 

In order to generate the code for a process, the user needs to prepare an input file, called process 
card in the following, which contains 

• process specific information, such as a list of initial and final state particles, their helicities 
(optional) and the order of the coupling constants; 

• scheme specific information and approximations, such as the regularisation and 
renormalisation schemes, the underlying model, masses and widths which are set to zero, 
the selection of subsets of diagrams, etc; 



Table 1. Part of the process card defining the initial and final state and the type and order of 
the coupling at LO and NLO 



process_path= qqgz 

in= u,u~ 

out= g, e-,e+ 

helicities= +- + -+,+ +,- + + -+,- + + 

order= QCD, 1,3 



• system specific information, such as paths to programs and libraries or compiler options; 

• optional information for optimisations which control the code generation. 

For example, to compute the QCD corrections for the process uu — )• gZ^ — )• ge~e^, the first 
part of the process card looks as shown in Table[TJ The first line defines the (relative) path to 
the directory where the process files will be generated. GoSam expects that this directory has 
already been created. Lines 2 and 3 define the initial and final state of the process in terms of 
field names, which are defined in the model file. Besides the field names one can also use PDG 
codes [39\ HO] instead. For more details we refer to |17j and the reference manual. 

3.3. Code generation 

If the process card is called gosam. in, it can be invoked by gosam.py: 

mkdir qqgz 
gosam.py gosam. in 
cd qqgz 

All further steps are controlled by the generated make files; in order to generate and compile all 
files relevant for the matrix element one needs to invoke 

make compile 

The generated code can be tested with the program matrix/test . f 90. The following sequence 
of commands will compile and run the test program: 

cd matrix 
make test . exe 
. / test . exe 

The last lines of the program output should look as follows. 



# 




LO: 


0.3450350717601E-06 


# NLO, 


finite 


part 


-10.77604823456547 


# NLO, 


single 


pole 


-19.98478948141949 


# NLO, 


double 


pole 


-5.666666665861926 


# IR, 


single 


pole 


-19.98478948439310 


# IR, 


double 


pole 


-5.666666666666666 



The printed numbers are, in this order, oq, co/oq, C-i/oq, c_2/ao and the pole parts calculated 
from the infrared insertion operator |41l H2] . One can also generate a pictorial representation of 
all generated diagrams using the command 

make doc 



which generates the file doc/process .ps 



3.4- Interfacing the code 

3.4- 1- Using the BLHA Interface The so-called Binoth Les Houches Accord (BLHA) [33] defines 
an interface for a standardized communication between one-loop programs (OLP) and Monte 
Carlo (MC) tools. GoSam can act as an OLP in the framework of the BLHA, such that the 
calculation of complete cross sections is straightforward. 

In general, the MC writes an order file, called for example olp_order . Ih, and invokes the 
script gosam.py as follows: 

gosam . py — olp olp_order . Ih 

The invocation of gosam.py generates a set of files which can be compiled with a generated 
make file. The BLHA routines are defined in the Fortran module olp_module but can also be 
accessed from C programs. 



3.4-2. Using External Model Files GoSam can also make use of model files generated by 
either Feynrules [43j in the UFO format [31j or by LanHEP [32j. The particles can be specified 
by their PDG code. Details about how to import thes model files are described in the Go- 
Sam reference manual. Precompiled MSSM_UFO and MSSM_LHEP files can also be found in the 
subdirectory examples/model. The examples directory also contains examples where UFO or 
LanHEP model files are imported. 



4. Examples 

The code generated by GoSam has been compared to a considerable number of processes 
available in the literature, as listed in Table [2j Many of these processes are also included 
as examples in the code, including reference values. 



W + 1 jet: transverse momentum of first jet 



- GoSam+Sherpa 

- MCFM 





(a) Transverse momentum of the leading jet for W +jet (b) Pseudorapidity of the leading jet for + jet production 
production at LHC with ^ = 1 TeV. at LHC with = 7 TeV. 



Figure 1. NLO calculation of W +jet production at LHC using GoSam interfaced with 
SHERPA via the Binoth Les Houches interface and compared to MCFM. 



As an example for te Binoth Les Houches Accord interface of GoSam we present results 
for the QCD corrections to W~ + 1 jet production, obtained by linking GoSam with 
SHERPA [59]. Furthermore, SHERPA offers the possibility to match NLO calculations with a 
parton shower |601 [6T] . Results for the transverse momentum and rapidity distribution of the 
leading jet are shown in Figs. 1(a) and 1(b) The comparison with MCFM [51 [B] shows perfect 
agreement. 
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Table 2. Processes for which GoSam has been compared to the literature. 



5. Conclusions 

We have presented the program GoSam which can produce code for the evaluation of one-loop 
matrix elements for multi-particle processes in an automated way. The program is publicly 
available at http : //projects . hepf orge . org/gosam/ and can be used to calculate one-loop 
amplitudes within QCD, electroweak theory, or other models which can be imported via an 
interface to LanHEP or FeynRules. Monte Carlo programs for the real radiation can be easily 
linked via the interface defined by the Binoth Les Houches Accord. 

The amplitudes are generated in terms of Feynman diagrams and can be reduced by unitarity 
based reduction at integrand level or traditional tensor reduction, or a combination of the two 
approaches. This feature makes the program extremely robust against numerical instabilities. 
The user can choose among different libraries for the master integrals, and the setup is such 
that other libraries can be linked easily. 

The calculation of the rational terms can proceed either together with the same numerical 
reduction as the rest of the amplitude, or before any reduction, using analytic information on 
the integrals which can potentially give rise to a rational part. Moreover, the GoSam generator 
can produce code for processes which include unstable particles, i.e. intermediate states with 
complex masses. 

GoSam is very well suited for the automated matching of Monte Carlo programs to NLO 
virtual amplitudes, and therefore can be used as a module to produce differential cross sections 
for multi-particle processes which can be compared directly to experiment. Thus GoSam can 
contribute to the goal of automating NLO corrections to an extent where NLO tools become a 
standard for data analysis at the LHC. 
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